Generating optical nonlinearity using trapped atoms 
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We describe a scheme for producing an optical nonlinearity using an interaction with one or more 
ancilla two-level atomic systems. The nonlinearity, which can be implemented using high efficiency 
fluorescence shelving measurements, together with general linear transformations is sufficient for 
simulating arbitrary Hamiltonian evolution on a Fock state qudit. We give two examples of the 
application of this nonlinearity, one for the creation of nonlinear phase shifts on optical fields as 
required in single photon quantum computation schemes, and the other for the preparation of optical 
Schrodinger cat states. 
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Our ability to perform quantum transformations on 
optical fields is hampered by the lack of materials with 
intrinsic optical non-linearities. While it is possible to 
circumvent this problem with schemes conditioned on 
photo-detection (see for instance 0, 0, EJ 0, H El)> 
ficiency problems are encountered both due to the in- 
herent nature of the schemes and efficiencies of current 
photo-detectors. In this paper we propose the use of the 
interaction with a simple atomic system, conditioned on 
high efficiency atomic measurements, to generate a near- 
deterministic (with probability approaching one) nonlin- 
earity on an optical field. This nonlinearity, together 
with linear transformations, is sufficient for generating 
arbitrary Hamiltonian evolution on a qudit formed from 
a truncated sequence of Fock states. 

The interaction between the single mode field and a 
two-level atom is described by the effective Hamiltonian 
H = K(a^cr~ + acr + ) where a\ a are the boson cre- 
ation and annihilation operators for the field mode and 
the operators cr + = \e)(g\, cr = \g)(e\ are the raising 
and lowering operators for the atomic state. The unitary 
transformation that acts when this interaction is applied 
for a time t is U(t) = exp[— ir(al<r _ + aer+)] where 
r = nt. 



If the atom is prepared in the ground state and found 
in the ground state after the interaction, the conditional 
state of the field is given by 




|2> = |e) 



FIG. 1: Level scheme for an effective two- level transition con- 
trolled by a stimulated Raman process. 
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On the other hand if the atom is prepared in the ex- 
cited state and found in the excited state after the in- 
teraction, the conditional state is given by T e (r)|V') = 
cos(tV aa^)\%l)). There is considerable practical advan- 
tage to using the excited state preparation rather than 
the ground state as there is always a signal for correct 
operation, however the analysis is the same, and we will 
concentrate on the operator 

We have in mind a quantum computing communica- 
tion protocol in which the optical field mode is derived 
from a transform limited pulsed field which is rapidly 
switched into the cavity mode containing the atomic sys- 
tems at fixed times determined by the pulse repetition 
rate. Similar systems have been proposed as a quantum 
memory for optical information processing Q. When 
the atomic measurement yields the required result the 
field may be switched out again for further analysis or 
subsequent processing through linear and conditional el- 
ements. 

Once the cavity field is prepared, we need to switch 
on the interaction with the atomic system. In order that 
we can switch this interaction at predetermined times we 
propose that an effective two level transition connected 
by a Raman process with one classical field and the quan- 
tised signal field, be used. A similar scheme has recently 
been proposed as the basis of a high efficiency photon 
counting measurement |3, li| • The process is also used 
in the EIT schemes for storing photonic information 0] 
and for quantum state transfer between distant cavities 
[lll |. The level diagram is shown in figure The nearly 
degenerate levels |1) and |2) are connected by a stim- 
ulated Raman transition to level |3). The detuning of 
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the Raman pulse from the excited state |3) is A, which 
is approximately the same as the detuning of the signal 
mode form the same transition. With these parameters 
the interaction strength is given by k = Qg/2A where f2 
is the Rabi frequency for the Raman pulse and g is the 
one photon Rabi frequency for the signal field. 

An advantage of using a stimulated Raman process of 
this kind is that the excited state |2) can be a meta- 
stable, long lived level. We thus do not need to consider 
spontaneous emission from this level back to the ground 
state. The readout of the atomic system may be achieved 
by using a cycling transition between the excited state 
1 2) and another probe level |4). Such measurements are 
routinely performed in ion trap studies |l2l | and can have 
efficiencies greater than 99%. 

In the following, we will prove that the ability to per- 
form the conditional transformation Q together with lin- 
ear transformations is universal on qudits. We will follow 
this with two examples. 

It is first instructive to consider a unitary operator of 
the form exp(i8VaJa). A nonlinearity of the form \/ at a, 
which we refer to as the square-root-number operator, 
turns out to be as good as a Kerr nonlinearity for uni- 
versal quantum computation in the infinite dimensional 
Hilbert space of a Harmonic oscillator . The concept 
of universal computation is isomorphic to the concept of 
universal simulation, which is the ability to simulate any 
arbitrary Hamiltonian evolution, to any degree of accu- 
racy, by combining the evolution due to a fixed class of 
Hamiltonian generators. 

To motivate this idea, consider what kind of nonlinear 
oscillator would correspond to the Hamiltonian H srn = 
\J a^a. The resulting Heisenberg equation of motion for 
the field amplitude operator is 



= at a - 



(2) 



which clearly indicates an oscillator in which the fre- 
quency is state dependent. The right hand side may be 
expanded as a power series in (at a) -1 / 2 to give 



d = -i(i(ata)-^-I (a ta)- 3 / 2 



+ . . . a (3) 



This corresponds to an oscillator for which the frequency 
of oscillation decreases with increasing energy. For com- 
parison, the Kerr nonlinear oscillator corresponds to an 
oscillator in which the frequency increases with increas- 
ing energy. It is thus clear that the square-root-number 
Hamiltonian will result in a rotational shearing of states 
in the phase plane of the oscillator in much the same way 
as occurs for the Kerr nonlinearity (l4|. 

We can go f urther by using the results of Braunstein 
and Lloyd [l3j. They considered the question of what 
Hamiltonians are universal for quantum computation in 
an infinite dimensional Hilbert space, such as that for a 
single mode of the radiation field. Their results show that 



Hamiltonians that are at most quadratic in the canonical 
momentum and position variables are not universal. For 
instance, in the case of a single mode field with annihila- 
tion and creation operators a, at, successive applications 
of Hamiltonians from the set of displacements, squeezing 
and rotations, 7Yii n = {za + z*a\ za 2 + z*(at) 2 , za)a} : 
can generate arbitrary linear canonical transformations 
in the phase space variables but no other transforma- 
tions. A universal set is easily obtained by adjoining al- 
most any Hamiltonian that is at least cubic in the canon- 
ical variables. A universal set of Hamiltonians could be 
made up from 7i\{ n together with the Kerr nonlinearity 
Hk = (at) 2 a 2 . Another choice is the cubic Hamiltonian 
H c = a) a(a + a') + he. We now show that the square- 
root-number operator VaJa together with the set 7in n 
can be used to simulate a cubic Hamiltonian, and thus 
is universal for quantum simulations in the Hilbert space 
of a single mode. 

Consider acting on the unitary operator generated 
by the square-root-number operator with a large dis- 
placement a: U(a,9) = D\a) exp(i6VaJa)D(a) where 
D{a) is the displacement operator given by D{a) — 
exp(aat — a* a). The resulting unitary operator may 
then be written as 



U(a, 



= e ie|Q| exp 



a* a 



a'ax(4>) 



2 2\a 



a 



8\a\ 



+ o(M- J ) 



(4) 



where x(<j>) = (ae -i * + he) with the phase of a. The 
first three terms in the argument of the exponential corre- 
spond to a second order Hamiltonian and simulate a dis- 
placement, a rotation, and a squeezing operation respec- 
tively. The fourth term however is a cubic term which 
is what we required for a universal set of Hamiltonians 
for a single mode. By using a linear Hamiltonian to mix 
several modes (for instance a beam-splitter), arbitrary 
multi-mode Hamiltonians can be constructed It is 

thus clear that we can use the square-root-number op- 
erator, together with an arbitrary linear transformations 
to perform universal computation. 

Now that we have shown the universality of the op- 
erator exp(idV ata) what can we say directly about 



cos 



at a ? The action of this conditional operator on 

the number state |n) is to multiply the state by the ampli- 
tude cos(0 v / n) = [exp(i8y/n) + exp(— i0\fn)\/2. Clearly, 
if we can restrict the interaction so that this amplitude is 
±1, then it is equivalent to the full unitary operator — and 
we might expect to exploit the nonlinear interaction in a 
similar way. Now, as it turns out, it is possible to 
choose 9 such that for a finite size computational space, 
\0)...\N), 



cos 



for n = 2 (2m + l) 2 
for n ^ 2(2m+ if 



(5) 
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Here m is a non-negative integer and we see a sign shift 
only on the states |2), |18), |50) ...\2(2m+ l) 2 ) .... We 
have chosen the first sign shift to occur on the |2) Fock 
state. This is an arbitrary choice and other terms can be 
shifted instead. 

For most computation schemes in which we are inter- 
ested, there is a finite number of basis states (there gener- 
ally always is a natural cutoff even in most CV schemes). 
In this case it is always possible to choose 9 to satisfy 
the above relations. This indicates that our operator 



cos 



is approximating an effective, highly non- 
linear unitary operator. For instance, if we consider only 
the qutrit subspace of |0), |1), and |2) then the operator 
acts equivalently to the Hamiltonian, H = which 
contains a Kerr nonlinearity. 

With this nonlinear Hamiltonian, it is possible in con- 
junction linear Hamiltonians to generate arbitrary multi- 
mode Hamiltonians. Hence we can in principle create 
any of the unitary operators required for universal qutrit 
computation. This argument extends to higher qudits in 
a direct manner. The potential problem however, is that 
the greater the number of states in the qudit space, the 
larger the value of 9 required to satisfy JHJ. 

It has recently been shown by Knill, Laflamme and 
Milburn (KLM) Q, that a conditional nonlinear sign 
shift (NS) gate on two photon states can be produced 
with passive linear optical elements and photo-detection. 
Such conditional nonlinear phase shifts can be used to 
perform two qubit operations for logical states encoded 
in photon number states. If such conditional gates are 
used to prepare entangled states for teleportation, ef- 
ficient quantum computation can be performed which 
with suitable error correcting codes can be made fault 
tolerant . Here we show that if the ancilla modes are 
replaced with a two level atom similar conditional non- 
linear phase shifts can be achieved by near determinis- 
tic post-selection on atomic measurements. The atomic 
measurements can be made with fluorescence shelving 
techniques which are very much more efficient than sin- 
gle photon counting measurements, thus reducing the 
need for new photon counting technologies inherent in the 
KLM scheme. In addition, the near-deterministic nature 
of our gate would also drastically simplify the implemen- 
tation of the KLM scheme. The ability to induce con- 
ditional nonlinear phase shifts on quantum optical fields 
also has applications to high precision measurements (see 
for instance 0). 

Our objective for the NS gate is to find a way to pro- 
duce the nonlinear phase shift defined by 




FIG. 2: A plot of the Q-function versus the two canonical 
phase space variables, for conditional state produced from an 
initial coherent state using, (a) a = 10, 9 = lOw and (b) 
a = 10, 9 = 5tt 



the state A o c o |0) + AiCi|l) + A 2 c 2 \2) where A = 1, 
A\ = cos(r) and A 2 — cos(v / 2t). Solutions for which 
(Ai,A 2 ) approaches arbitrarily close to (1,-1) can eas- 
ily found, for examples see table [I] It is also possible 



r 


A Ai A 2 


6.5064 
37.73742 
219.918 


1 0.97519 -0.97516 
1 0.9992663 -0.9992665 
1 0.999979 -0.999978 



TABLE I: High-probability results for a single atom initially 
in a ground state and measured in a ground state after an 
interaction time r. 

to employ several atoms, each addressed independently, 
and find solutions of high probability using the same tech- 
niques. For example with two atoms, one initially pre- 
pared in the ground state and the second prepared in 
the excited state, the conditional state given that both 
atoms are found in their initial state after the interac- 
tion is \4>ge) — Y g(T~i)~£ e (T 2 )\i(}} and again values of the 
interaction times t% and t 2 can be found which again 
perform the NS gate with high probability. For instance 
with n = 37.79300921 and t 2 = 197.78109842, then 
I Ai| = |A 2 | = | A3 1 = 0.990321935 and the required phase 
shift is performed on |2). 

In any real experiment, the desired interaction time 
can be calibrated by placing the phase shift in one arm 
of a Mach-Zehnder interferometer, with two single pho- 
ton inputs, and examining the interference fringes as a 
function of the interaction time. 

Now let us consider the case of the field in an initial 
coherent state \a). We will show that the conditional 
transformation, cos(9V ata), generates a coherent super- 
position of coherent states, a so called Schrddinger cat 
state. A convenient phase space representation of the 
conditional state is the Q-function defined by 



Q(P) = \(f3\cos(9VaUi)\a)\ 2 



co\0) + Cl \l)+ c 2 \2) -» c |0> + Cl |l) - c 2 |2) (6) 

where the Cj are complex amplitudes satisfying X) l c il 2 = 
1. This is in the form of the conditions in equation 
so consider applying the measurement operator T g , on 
the general two photon state. This will leave us with We thus first consider the amplitude function A(f3) 
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(f3\ exp(i6>V aJa)\a). Based on the semi-classical expec- 
tation that this unitary transformation describes an oscil- 
lator with an energy dependent frequency, we anticipate 
that we need only consider the Q-function on the curve 
|/3| = \a\. With this in mind we put (3 = \a\e % ^ and the 
Q-function amplitude is then given by 



y i(8s/n—(j)n) 



(7) 



n=0 



where p n (a) — e~' Q ' 2 ^-. We now assume that \a\ >> 1 
and approximate the Poisson distribution, p n (a) with a 
Gaussian 



Pn (a) « (27r|a|)- 1/2 exp 



(n-|a|) 



2 1 



2 a 



(8) 



We can then replace the sum with an integral over the 
variable y = n — \a\. Under the assumption |a| >> 1 
we find the integrand can be approximated as a general 
Gaussian and thus the integral is given by 



A{p) = e"^ |Q| exp 



i(tf-2|a|0) 2 



(9) 



Clearly this distribution is peaked at (f> = 9/(2\a\). If we 
choose the interaction time so that 9 = \a\ir, we expect 
the state to be localised on the positive imaginary axis 
in the phase plane of the Q-function. Similarly for the 
same parameters (f3\ exp(i6'V aJa)\a) will be localised on 
the negative imaginary axis in the phase plane of the 
Q-function. It then follows that for the full conditional 
operator, cos(6'v / aJa), the state has two components lo- 
calised symmetrically about the origin on the imaginary 
axis. In figure we show the Q function as a function 
of the two canonical phase space variables x and p for 
several choices of a and 9. 

We observe that if we load a coherent state into the 
cavity, and repeat the conditioning measurement proce- 
dure outlined previously, the conditional state of the field 
will be prepared in a state which is close to a cat state. 
Note however that these cats are not parity eigenstates 
as the conditional interaction cannot change the photon 
number distribution. Similar cat states are produced by 
a Kerr nonlinearity 0, 0] . 

We now estimate some typical values for the param- 
eters. In a recent experiment a similar stimulated Ra- 
man process was observed using single rubidium atoms 
falling through a high finesse optical cavity ^3] • The fol- 
lowing parameters are typical of that experiment: g = 
2tt x 4.5 MHz, Vt = 2tt x 30 MHz and A = 2tt x 6 MHz. 
This gives a coupling constant of the order of 70 Mhz. 
To achieve effective interaction constants of the order of 
those in the table [i] requires interaction times of the or- 
der of 0.1 — 5 /is. In this paper we have neglected cavity 
decay which obviously needs to be kept small over similar 
time scale, which while difficult is not impossible. 



We have shown how the ability to do very efficient 
measurements on single atoms trapped in an optical cav- 
ity can be used to implement nonlinear conditional phase 
shifts on the intra-cavity field. By carefully choosing the 
interaction time, a nonlinear interaction can be imple- 
mented with near unit probability, that, together with 
linear transformations, is universal for simulating an in- 
teraction on qudits. For small qudits the required inter- 
action time can easily be found numerically. If the field 
state can be carefully switched in and out of the cavity, 
the method can be used to implement near-deterministic 
nonlinear gates for quantum optical computing. For in- 
stance, the scheme can be used to implement a nonlinear 
sign shift gate, which thus provides a path to quantum 
computation with logical qubits encoded in photon num- 
ber states. It can also be used to conditionally generate 
coherent superpositions of coherent states and thus can 
provide the key resource for quantum computing with 
coherent states 
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